Volume 102, Number 6, November-December 1997 

Journal of Research of the National Institute of Standards and Technology 

[J. Res. Natl. Inst. Stand. Technol. 102, 677 (1997)] 

Current Distributions in Quantum Hall 

Effect Devices 



Volume 102 



Number 6 



November-December 1997 



M. E. Cage 

National Institute of Standards and 

Technology, 

Gaithersburg, MD 20899-0001 



This paper addresses the question of how 
current is distributed within quantum 
Hall effect devices. Three types of flow 
patterns most often mentioned in the lit- 
erature are considered. They are: (1) skip- 
ping orbits along the device periphery 
(which arise from elastic collisions off 
hard-walled potentials); (2) narrow con- 
ducting channels along the device sides 
(which are presumed to be generated 
from confining potentials); and (3) currents 
distributed throughout the device (which 
are assumed to arise from a combination 
of confining and charge- redistribution 



potentials). The major conclusions are that 
skipping orbits do not occur in quantum 
Hall effect devices, and that nearly all of 
the externally applied current is located 
within the device interior rather than along 
the device edges. 
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1. Introduction 



There is considerable discussion and disagreement in 
the literature about how the current is distributed within 
quantum Hall effect devices. Some authors assume that 
the current is confined to skipping orbits along the 
device periphery, others believe it is confined to narrow 
edge- state channels along the device sides, and still oth- 
ers believe it is distributed throughout the device inte- 
rior. 

After a brief discussion of the quantum Hall effect, 
this paper investigates the nature of current patterns 
arising from three models that cover the above possibil- 
ities: (1) skipping orbits; (2) narrow conducting chan- 
nels due to confining potentials along the device sides; 
and (3) distributed currents due to either a charge-redis- 
tribution potential within the device interior or to a 
combination of a charge-redistribution potential and 
confining potentials along the device sides. It is pointed 



out when the predictions of the three models agree or 
disagree with experimental results. 



2. Integer Quantum Hall Effect 

The integer quantum Hall effect [1-3] requires a fully 
quantized two-dimensional electron gas (2DEG). On an 
even integer Hall plateau, conducting electrons of the 
2DEG have completely filled all the allowed spin-down 
and spin-up states of the lowest Landau levels, and none 
of the next Landau level. Negligible dissipation occurs 
within the interior of the 2DEG in the Hall plateau 
regions of high-quality devices operated at low enough 
currents. Within these regions the quantized Hall resis- 
tance Rh of the i\h plateau has the value /?h(0 = h/(e^i), 
where h is the Planck constant, e is the elementary 
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charge, and / is an integer. In GaAs the quantum number 
for the first filled Landau level is / = 2. 



3. Skipping Orbit Currents 

Electrons of the 2DEG execute cycloidal motion 
when in the presence of an applied perpendicular mag- 
netic flux density B . The cyclotron radius is Tc = (^/ 
eBY^ for orbits of the first Landau level, in which case 
each electron of the 2DEG has trapped a magnetic flux 
quantum hie . If: (a) the device is homogenous; (b) there 
is no applied current /so between the source S and the 
drain D of the device; (c) the magnetic flux density B is 
uniform; and (d) the device boundaries are represented 
by hard-walled confining potentials, then semiclassi- 
cally the electrons are uniformly distributed throughout 
the device interior. The cyclotron orbital velocities vec- 
torially tend to cancel everywhere within the device, 
except near the device periphery where there are skip- 
ping orbits due to elastic scattering from the hard- 
walled potential, as indicated in Fig. 1. An electron 
current therefore circulates around the device boundary 
in the absence of an applied current, thereby generating 
a measurable magnetization [4]. 

It is straightforward to estimate the magnitude of this 
skipping orbit current for a homogenous device when 
/sD = 0. For a GaAs/AlGaAs heterostructure device, in 



which the / = 2 plateau occurs at a typical magnetic flux 
density B = 12.3 T, the cyclotron radius Tc = {/ileBy^ is 
7.3 nm; the orbital angular frequency (jl>c = eBIm^ is 
3.2 X 10'^ /s, where m* is the reduced mass of the 
electron (0.068 times the free electron mass m^\ and the 
orbital velocity Vc = Wc'^c is 2.3 X 10^ m/s. The electrons 
travel a distance Ivc along the device boundary in a time 
t which is one-half the cyclotron orbital period Tc, so 






2/c 



1277 

■2c^. ' 



(1) 



where the orbital frequency /c is 5.1 X 10^^ Hz and t is 
9.8 X 10 '"^ s. The average electron velocity for skip- 
ping orbits along the device periphery is thus Vso = 2rjt, 
or 1.5 X 10^ m/s. The average density of the 2DEG is 
n^ = ieB/h, or 5.9 X 10'^ /m^; and the skipping orbit 
current density is Jso = n^evso, or a very large 1 42 A/m. 
The skipping orbit current is therefore 



ie^B 



(2) 



which is similar to the expression 1^^ = ie^Blilum^) 
stated without proof by Thouless [5]. Thus, according to 
Eq. (2), 1.0 |ulA of skipping orbit current circulates 
around the device periphery for the / = 2 plateau at 12.3 
T in the absence of an applied current. 




V 



so 



JSD=0 



•so 



Fig. 1. Skipping orbits of tlie conducting electrons around part of tlie device peripliery wlien the 
magnetic flux density B points into the figure. The average velocity of the skipping orbits along the 
device boundary is Vso, and the skipping orbit current is Uo. There is no externally applied current 
/sD- The lightly shaded region represents a mesa etch down below the 2DEG, while the darker 
shaded region is an ohmic contact to the 2DEG. 
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What happens if there is an externally applied cur- 
rent? Consider the case where all the applied current 
executes skipping orbits. Good quantization of the / = 2 
quantum Hall voltage Vh, and small values of the longi- 
tudinal dissipative voltage Vx, have been observed for 
applied currents as large as 200 fxA [6, 7]. The above 
assumption would yield an enormous skipping orbit cur- 
rent density Jso = hJrc = 2.7 X 10"^ A/m, and an average 
electron velocity Vso = J^J{n^e) = 2.9 X 10^ m/s, that is 
9.7 % the speed of light in vacuum and 35. 1 % the speed 
of light in GaAs. This current density and velocity is 
unrealistic, but it is sometimes stated [5, 8-10] that for 
small currents, where Vh is less than the Landau level 
spacing ^coje, that the applied current is an edge cur- 
rent, and that this edge current is concentrated within a 
cyclotron radius of the device boundary, and therefore 
undergoes skipping orbits. These authors assume that 
for larger currents, where Vh > ^(oje, part of the ap- 
plied current is a skipping orbit current and part is a bulk 
current within the device interior. 

Figure 2 shows the skipping orbit flow patterns for 
the case of a magnetic flux density pointing into the 
figure, and applied currents ± /sd small enough to sat- 
isfy the assumption that the current is contained entirely 



within skipping orbits. The thin lines represent an inter- 
nally induced skipping orbit current, as in Fig. 1 , while 
the thick lines are for an internally induced current plus 
the externally applied current /sd- The applied current is 
only along one side of the device in Fig. 2 because: (a) 
the Lorentz force ev^B^ is equal and opposite to the 
Coulomb repulsive force eEy everywhere within the 
2DEG; (b) Ey= — Wy of the hard-walled confining po- 
tential has the opposite sign on each side of the device; 
and (c) only one side of the device has the appropriate 
sign of Ey for the applied current direction. 

There are two serious problems with the flow patterns 
of Fig. 2. The first problem is that current enters and 
exits the bottom corners of the source and drain contacts 
for one current direction and the top corners for the 
opposite current direction, whereas there is clear exper- 
imental evidence [11-14] that the current enters and 
leaves opposite corners of the device, and that these 
corners remain the same when the current is reversed. 
The second problem is that the current enters and exits 
every potential contact in the figure. (Even if the trans- 
mission coefficient is not unity on the first attempt to 
enter the contact, the reflected electrons skip along the 
contact and eventually enter it.) However, there is exper- 
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Fig. 2. Current-carrying paths through the device if all the current were carried via skipping orbits. 
The thin lines are for an internally induced skipping orbit current as in Fig. 1. This current exists even 
when /sd = 0. The thick lines represent internally induced current plus an externally applied current 
/sd- Figure 2(a) is for an external current of electrons entering the source contact S and exiting the 
drain contact D; Fig. 2(b) is for the opposite current direction. The magnetic flux density points into 
the figure in the positive z direction. 
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imental evidence [15, 16] that the resistivity of the 
2DEG is less than 5 X 10"^ fl at 1.2 K and less than 
10" '° fl at 0.4 K. The latter value corresponds to a 
three-dimensional resistivity of less than 10 ^^ fl cm, 
which is much smaller than the resistivity of any other 
non- superconducting material. The contacts themselves, 
however, are resistive in high magnetic fields. Thus, 
under steady-state dc conditions, the current will take 
the path of least resistance and avoid the potential con- 
tacts. (This assumes that no current is drawn through the 
potential probe contacts by an external measurement 
system, and that the time is long enough after current 
reversal for the device potentials to have reached equi- 
librium.) 

Skipping orbits do not seem feasible in the quantum 
Hall effect for the two reasons just given, and for rea- 
sons that will be presented in Sees. 4.3 and 5.3. 



4. Narrow Conducting Channel Currents 

Even if skipping orbits do not occur in the quantum 
Hall effect, current will still circulate around the device 
periphery because a confining potential must exist to 
prevent electrons of the 2DEG from spilling out of 
mesa-etched devices when the etching extends below 
the depth of the 2DEG. We will see in this section that 
confining potentials along the sides of the device create 
narrow channels, along which, current can flow. 



4.1 Eigenstates of the 2DEG 

Reference [17] explains how the confining potential 
arises from electron surface charges on the side of the 
mesa and a charge-depletion region near the device 
edge. Note that this potential is negative. Therefore we 
assume a negative confining potential Vc of finite spatial 
extent A on each side of the device. The energy of an 
electron in this confining potential is then 
Sc = qVc = — eVc. Figure 3 is a schematic drawing of 
the electron energies, plotted as a function of j, across 
the device width w for a confining potential of greatly 
exaggerated spatial extent A . (For convenience, the fig- 
ure shows a linearly-shaped confining potential, but we 
will use a more realistically shaped parabolic confining 
potential in the calculations.) No current exists within 
the device interior in Fig. 3 because the potential gradi- 
ent is zero between —A and A. 

The electrons in the 2DEG occupy unique quantum 
eigenstates, indicated as circles in Fig. 3. The usual 
method of defining these states is to represent their 
wavefunctions in the Landau gauge as normalized prod- 
ucts of Hermite polynomials across the device multi- 
plied by plane waves propagating down a length Ljc of 
the device [ 1 8-20] . Let us consider only even-/ quantum 
Hall plateaus. The energy eigenvalue Sn of each state in 
Landau level A^ is then 



S^iyo) = (A^+ 1 )^Wc + eyoEiyo) + ^ m^v.'Cjo), (3) 
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Fig. 3. Schematic drawing of tlie energies of allowed eigenstates across the device width w for 
a linearly shaped confining potential of greatly exaggerated spatial extent A and no applied current 
/sD- Only a few eigenstates are indicated. The jc -axis points along the center line of the device. The 
electrochemical potential jjl is the same on both sides of the device for this case. Eigenstates of 
the lowest Landau level 1 are filled between jmax = ~ Jmin- Shaded circles are for occupied states 
that do not contribute to the current, black circles are current-carrying occupied states, and open 
circles are unoccupied states. 
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where 



Jo = 



■■{vJoj^+ £IK) 



(4) 



is the center-of-mass position of each state undergoing 
cycloidal motion, — wl2 < jo < wl2\ Vxiy) = E{y)IB is 
the electron drift velocity down the device; £b = {M 
cBY^ is the magnetic length, and is equal to the cy- 
clotron radius r^ for the first Landau level; and 
k^ = lirNk/Ljc is the wavevector for the state located at 
position Jo with an associated positive or negative in- 
teger quantum number Nk. The eigenstates are repre- 
sented by the quantum numbers (N,Nk), and the wave- 
function for each state is 



^N,N,(x^y)-- 



1 



(LxV' 



^ilnNj^x/Lj^ 



1 



(2*A^!)'' 



X 



(^4^)' 



.e-(.-.o)^/24H^[(j,_j,„)/4], 



(5) 



where H^Ky — yoV^s] is a Hermite polynomial. 

The eigenstates, represented as circles located at cy- 
clotron center-of-mass positions jo, are shown in Fig. 3 
for the first (A^= 0), (i = 2) Landau level and the second 
(A^= 1), (i = 4) Landau level. Only a few of the eigen- 
states are indicated . The sp atial extent of each eigenstate 
[18-20] is ± €bV2N-\- 1, which is equal to 2rc for the 
first Landau level; the spatial separation between adja- 
cent states in a constant electric field is Ajo = 2i:€llLj, 
from Eq. (4); and the energy separation between adja- 
cent Landau levels is ^Wc- 

The magnetic flux density has been adjusted in Fig. 
3 so that all the allowed states of the first Landau level 
are filled and none of the second. The Landau level is 
midway across a mobility gap [21]. Only localized states 
due to imperfections and impurities are being filled. 
These localized states (not shown in the figure) do not 
affect the quantum Hall voltage because the states are 
stationary [22-23] (non-conducting). Since there is no 
applied current, and therefore no Hall voltage, the Fermi 
energy S^ is constant across the device width, and is 
located halfway between Landau levels, as indicated by 
the dotted line in Fig. 3. Under these conditions, states 
of the lowest Landau level are occupied up to the Fermi 
energy S^ = fioy^, and no states are occupied in the sec- 
ond Landau level. The occupied states are located be- 
tween jmax, and j^in, and in this case | j^ax I = I Jmin I ■ 
Shaded circles are for occupied states that do not con- 
tribute to the current because Ey is zero; black circles 
represent current-carrying occupied states where 
Ey 4" 0; and open circles are unoccupied states. 

Voltage probes located along the sides of the device 
measure the electrochemical potential ^t of the 2DEG 



[24]. We assume that the probe potential is that of the 
nearest occupied conducting state on that side of the 
device, i.e., the state at j^ax or jmin- We also assume that 
the conducting states at j^ax or jmin are far enough from 
the mesa-etched sides that the shape of the confining 
potential is not significantly altered by the boundary 
conditions, and that the potential and the electric field 
are continuous across the mesa-etched interface. In the 
special case of Fig. 3 where /sd = 0, ^/Ar = e/^L = ^f on 
either side of the device. 

4,2 Confining Potential 

A current circulates around the device periphery 
within the narrow regions containing the black-circle 
eigenstates of Fig. 3 when there is an external magnetic 
field. The existence of this current is predicted by edge- 
channel models [8, 24—31], and is verified by experi- 
ment [32-37]. (Edge channel is a better description of 
this phenomenon than edge state because each conduct- 
ing channel is composed of many eigenstates.) Most 
edge-channel models recognize that the confining po- 
tentials have a finite extent, but then assume that the 
confining potentials are hard-walled at ± w/2 when 
making calculations, and that j^ax = >v/2 — Tc = — jmin- 

We consider the more realistic case where the confin- 
ing potentials have a finite spatial extent, rather than 
hard walls. For simplicity, the confining potentials 
shown in the schematic drawing of Fig. 3 are linear, but 
we will assume parabolically-shaped confining poten- 
tials in the calculations. These parabolically-shaped 
confining potentials arise from homogeneous charge- 
depletion regions of spatial extent A, The confining 
potentials of Fig. 3 have origins at j = A = w/2 — A and 
j= — A= — w/2 + Zi. The value of A has been deter- 
mined in an experiment by Choi, Tsui, and Alavi [38]. 
They used one-dimensional localization theory to evalu- 
ate conduction in narrowly-constricted channels of 
GaAs/AlGaAs heterostructures at very small magnetic 
flux densities, and found that A = 0.5 |ULm ± 0.2 juim. 
(This result is consistent with observations that devices 
having widths less than 0.5 |ULm cease to conduct alto- 
gether at temperatures below 4.2 K [38], even at high 
magnetic flux densities [39] where the 2DEG penetrates 
into the depletion region.) We choose the value of the 
charge-depletion depth to be their average value 0.5 |ULm. 
The remaining parameter to define the confining poten- 
tial is Vm, the value of the confining potential at ± w/2. 
The confining potential exists even at zero applied cur- 
rent and zero magnetic flux density. Also, there are 
many impurity states between the valence and conduc- 
tion bands. A reasonable value of Vm is one-half of the 
separation voltage between the valence and conduction 
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bands, which is V^ = 0.75 V in GaAs at 1 K [40]. This 
value is comparable to the 0.8 V used by Choi, Tsui, and 
Alavi [38]. 

The equations for the parabolic confining potential Vc 
and its electric field Ec= — VVc are therefore 

Vc(y)=-a(y- \f and E,(y) = 2a(y - \) (6a) 



for A ^ J ^ ^ , 



(6b) 



Vc(y) = and E,(y) = 
for - A < J < A , 

V,(y)=-a(y-^\f and E,(y) = 2a(y -^ \) (6c) 



for -^< -};< - A, 



where a = VJA"- = 3.0 X 10" V/m' for A = 0.5 |ULm and 
Vm = 0.75 V, and 



A = f-., (7, 



4.3 Confining Potential at ho = |xA 

Given the above equations and the values of A and Vm, 
much information can be deduced about the current 
distribution when /sd = |ulA and the magnetic flux 
density is adjusted to be halfway between Landau levels 
and on the i = 2 quantum Hall plateau, as is the case in 
Fig. 3. Under these conditions 



^cCymax) = -^ = - eVciym^) = ^« (jmax " A )^ 



■ e J2 ( Jmax - 2 + ^ 



(8) 



if we measure the energy ^^cCVmax) relative to 6^c(0). 
Thus 



and 



Jmax = - Jmin = 199.559 fXm 



= 0.441 [xm 



(9) 



(10) 



for a 400 |ULm wide quantized Hall resistance standard 
device at a magnetic flux density of 12.3 T. The electric 
field E,(y) from Eq. (6a) is 3.6 X 10^ V/m at 



Jmax = 199.559 |ULm. The current density along thex-axis 
at position y is Jdy) = o-^E^iy) + o-^Edy) = o-^Edy) in 
the absence of significant dissipative scattering [41], 
where the off-diagonal conductivity tensor component 
is a^ = ie^lh=\IRn=\l\2 906.4 VL for the i = 2 
plateau. /cCVmax) = 27.5 A/m for this case. 

The total current carried by the occupied states of the 
right-hand side (rhs) confining potential is 



/c(rhs)= J^{y)dy= (T^E,{y)dy = 



Rn 



[K(jmax)-K(A)]: 



Rn ' 



(11) 



where KCymax) = ~ ^(jmax ~ A)^. Similarly for the left- 
hand side (Ihs), 

-A -A 

/c(lhs)= J /c(y)dj= J o-^c(y)dj = 






(12) 



where KCymm) = — ^ (jinm + A )^ It follows from Eqs. 
(8), (11), and (12) that 

for the 12 906.4 fl, / = 2 plateau at 12.3 T and /sd = 
|ulA. This rather large 0.8 fxA edge-channel current cir- 
culating around the device is comparable to the 1.0 |ulA 
current obtained in Sec. 3 for skipping orbits. Note from 
Eq. (13) that the current is independent of the device 
width w and of the confining potential parameters A and 
Vrn if w > 2A and /sd = |ulA, whereas j^ax and y^m 
depends on w, zl, Vm, and /sd- 

We see from Eq. (10) that the maximum spatial extent 
of the current-carrying states is 60 times farther away 
from the sides of the device than predicted for the 7.3 
nm radius skipping orbits; thus there is no need to in- 
voke skipping orbits. Indeed, if there were skipping or- 
bits, the electric field Efy) calculated from Eq. (6a) 
would be an enormous 2.96 X 10^ V/m at jmax = w/ 
2 — r^. Also, the current density JJ^I2 — r^ would be 
a very large 229 A/m, even for this case with no applied 
current. As a final argument against skipping orbits, 
note in Fig. 3 that there are unoccupied Landau eigen- 
states (and localized states not shown in the figure) at 
j^ax and j^i^. If the states at j^^x and y^.^ executed skip- 
ping orbits, then scattering into unoccupied states could 
occur at every reflection, and this scattering need not be 
elastic. Hence dissipation could occur. The electron 
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current paths will surely adjust to minimize this dissipa- 
tion. 

Next, consider the case in Fig. 4 where /sd = |ulA 
and two Landau levels are occupied, i.e., the case for the 
i = 4 plateau. The energies of the highest-filled eigen- 
states of the first and second Landau levels are (3/2)^a)c 
and (l/2)^a)c, respectively. We can use Eq. (8) to calcu- 
late Jmaxl = - Jminl aud Jrnax2 = " Jmin2 for £L typlcal 5.5 T 

magnetic flux density on the i = 4 plateau. The values 
are 199.568 [xm and 199.540 juim. The electric fields 
and current densities at the two values of j^ax are 
4.1 X 10^ V/m and 63.6 A/m, and 2.4 X 10^ V/m and 
36.7 A/m, respectively. Three-fourths of the 1.45 |ulA 
current circulating around the device periphery with no 
applied current is due to electrons in the lowest Landau 
level. 

We can also estimate the spatial separation of the two 
conducting edge-channels for the first two Landau lev- 
els: they are separated at the Fermi level energy by 
Jmaxi ~ Jmax2, aud thc total Spatial extents of th eir wave - 
functions are approximately €bC^2N-\- 1 + V2A/''+ 1). 
Their separations for this example (where N=0, N' = 1 , 
and B = 5.5 T) are 8 % farther apart than the spatial 
extents; thus we predict that the edge-channels are phys- 
ically separated from each other. Buttiker [28] used a 
different approach in his Eq. (56) to arrive at a similar 
result. 



4,4 Confining Potential at /sd ^ fJiA 

Figure 5 shows the situation for an applied current /sd, 
assuming all the current is within the confining potential 
regions A to j^ax and — A to j^in- (This is a reasonable 
assumption for very small applied currents in the pi- 
coampere and nanoampere ranges.) The value of j^ax 
increases with increasing /sd, and additional eigenstates 
are occupied on the right-hand side of the device. Fewer 
states are occupied on the left side, and | jmin I de- 
creases. The chemical potentials are now different on 
the two sides of the device, and the Hall voltage is 
Vh = ^h/sd = (Ml — Mr)- If the confining potentials 
were linear we could determine the values of j^ax and 
jmin for any small applied current because one-half the 
Hall voltage would appear on each side of the device. 
However, we require more information for other poten- 
tial shapes because a range of j^ax values exists for 
which corresponding values of jmin can be obtained that 
also provides the correct quantum Hall voltage. 

Figure 6 is a schematic of the current-carrying paths 
through the device if all the current is carried via edge- 
channels generated by confining potentials, as is the 
case in Fig. 5. The paths are for ± /sd and ± B . Current 
circulates around the device, and it would be equal and 
opposite on either side of the device if /sd = |ulA, as in 
Fig. 3. More of the current is carried by the thick-line 
paths when there is an applied current. The applied 
current /sd enters and exits opposite corners of the 
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Fig. 4. Similar to Fig. 3, but two Landau levels are filled and there are now two values of jmax and two values 

of Jmin = - Jmax- 
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Fig. 5. Similar to Fig. 3, but with a small value of applied current /sd- The chemical potential is now 
different on the two sides of the device. 



device, in agreement with experiment [11-14]. These 
opposite corners are interchanged on magnetic field re- 
versal. There are no skipping orbits in Fig. 6 because 

Jmax < (w/2 - Tc). 

I suggest in Fig. 6 that the dc current takes the path of 
least resistance and avoids the potential contacts once 
the device reaches steady-state conditions because the 
resistivity of the 2DEG is so much smaller than the 
resistivity of the potential contacts [15, 16]. Perhaps, 
however, there is a physical requirement for the current 
to enter the potential contacts, as assumed in the 
Landauer-Buttiker formalism [26-28] with transmission 
and reflection coefficients at the contacts. 

We assume in Fig. 5 that Vh « ^(oje and /sd « ^(oj 
eRii since there is no significant electric field or current 
in the interior between — A and A . For the i = 2 plateau 
at 12.3 T this means that /sd « 16 |ulA. This condition 
is easily satisfied in edge-channel experiments [32-37] 
where the current is typically less than 50 nA. This 
condition is also consistent with the experiment of Kane, 
Tsui, and Weimann [42] where they observed a change 
in the behavior of the longitudinal resistance Rj, between 
filamentary and bulk-like currents at 1.5 |ulA: at 1.5 |ulA, 
Vh = ^(oje in their device, and Rj, then scaled inversely 
with the device width w, as predicted from the bulk 
current condition Rjc = pjocLjc/w, where pjoc is the resistiv- 
ity and Ljc is the length down the device between V^ 
probes. 

This section has dealt with zero or very small applied 
currents, such that eVn « ^(o^. That is often not the 
case. Indeed, we have used applied currents as large as 
200 fxA and still observed reasonably good quantization 
of the Hall voltage Vh and small values of the longitudi- 
nal dissipative voltage V^ = RxIsd [6, 7]. A 200 |ulA cur- 



rent for the i = 2 plateau yields eVn = 123 ^(Oc, which 
means that the assumption of Fig. 5 that current does not 
exist in the device interior no longer holds. 



5. Distributed Currents 

There is experimental evidence [13, 42-A6] that sig- 
nificant current exists within the device interior when 
^Vh > ^(Oc. This is often referred to as a bulk current. 
Some of the experiments [43^6] used contacts within 
the device interior, leading to concern that the contacts 
perturbed the current distribution, and speculation that 
the apparent bulk current was really due to edge-chan- 
nels at each internal contact. However, part of the exper- 
iment of Kane, Tsui, and Weimann [42] used no internal 
contacts, and the experiment of Fontein et al. [13] was 
a contactless measurement. Those two experiments 
clearly indicate the existence of internal currents. This 
section will therefore consider the case of current within 
the device interior. 

5.1 Charge-Redistribution Potential 

Initially we will ignore the confining potentials, but 
they will then be included in Sec 5.2. Several theories 
have been used for bulk currents, such as a classical 
electrodynamics model using local conductivity tensors 
[47, 48], and percolation down the device along paths of 
constant potential [49-52] . A model will be used here in 
which the applied current /sd induces a potential distri- 
bution within the device that extends across most of the 
device width w . 
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Fig. 6. Current- carrying paths through the device if all the current is via edge-channels generated by the 
confining potentials, as in Fig. 5. The thick lines represent the side of the device where more of the current 
flows. There is less current through the thin lines. The magnetic flux density points into the device in the 
positive z direction in Figs. (6a) and (6b), and out of the device in Figs. (6c) and (6d). 



The Lorentz force exerted on the conducting elec- 
trons of the 2DEG causes an increase in the density of 
electrons on one side of the device and a decrease on the 
other side. Thus, there are deviations, — e8a(y), from 
the average surface charge density — ens= — ie^Blh of 
the 2DEG across the device width. The charge-redistri- 



bution — edaiy) can be represented as a sequence of 
line charges, where the sequence is across the device in 
the ± y directions and the line charges point along the 
device in the ± x directions. A logarithmic charge-redis- 
tribution potential Vriy) across the device results from 
this sequence of line charges. 
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MacDonald, Rice, and Brinkman [53] expressed 
Vr(y) self-consistently in terms of the charge-redistribu- 
tion. Riess [54] extended this potential to a 2DEG with 
finite thickness. Thouless [55] found an analytic loga- 
rithmic approximation of Vr(y) far from the sides of the 
device. The charge-redistribution potential is infinite at 
the physical device edge, so Beenakker and van Houten 
[29] approximated the near-edge behavior by introduc- 
ing a cut-off near the device side and a linear extrapola- 
tion to the edge. Balaban, Meirav, and Shtrikman [56] 
used a quadratic extrapolation near the device sides, and 
a cut-off that was the same at both edges of the device 
and at all currents. Their cut-off distance from the 
device edge was the magnetic length €b = (^/eBY^^ 
(which is the cyclotron radius Tc for skipping orbits of 
electrons in the first Landau orbit). Cage and Lavine 
[17] used the same form for the potential as Balaban et 
al. [56], but a different geometrical factor and the very 
different cut-off values 5max = w/2 — y^ax and 8mm = w/ 
2 + jmin (which differ on either side of the device, vary 
with applied current, and depend on the magnetic flux 
density direction). No extrapolation to the device edges 
was used; the occupied Landau eigenstates were as- 
sumed to be far enough from the device edges to be 
unaffected by the conditions that the potential and elec- 
tric fields are continuous across the boundary at the 
mesa edge. 

The charge-redistribution potential of Cage and 
Lavine [17] is 



, (14) 



Vr(y) = 


IrRn 

2 


^^ Jmax + w/2 
_ w/2 - Jmax. 


-1 

In 


y + w/2 
y — w/2 




- 


2 < Jmin ^y^ 


Jmax 


w 
<2 



for 



where 



/r = /sD-/c(rhs)-/c(lhs). 



(15) 



/c(rhs) and /c(lhs) are defined by Eqs. (11) and (12), and 
are zero if there is no confining potential. The geometry 
factor G is 



G(w,ymax) = 



In 



Jmax + w/2 



w/2 - ym 
The charge-redistribution electric field E^ = — VK is 



(16) 






■y 



(17) 



Figure 7 is a schematic drawing of the energy 
Sr = qVr = — eVr across the device width w for two 
Landau levels when using the charge-redistribution po- 
tential defined by Eqs. (14) and (15). A few energy 
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Fig. 7. Schematic drawing of the energies of allowed eigenstates of 
two Landau levels for the applied current-induced charge- redistribu- 
tion potential defined by Eqs. (14) and (15). Only a few eigenstates are 
indicated. Eigenstates of the lowest Landau level are filled between 
ymin and y^ax- Black circles are current-carrying occupied states that 
contribute to the chemical potentials ^r and /jll on either side of the 
device. Open circles are unoccupied states. 

eigenstates are also shown. Vr(y) is infinite at ± w/2, 
but that does not matter because there are no occupied 
eigenstates beyond jmax and ymm- Only the occupied con- 
ducting states of the first Landau level between y^ax and 
ymin (indicated in black) contribute to the chemical po- 
tentials on either side of the device and to the quantum 
Hall voltage. The potential is therefore finite and well- 
behaved in the region of interest. 

This defines a realistic charge-redistribution poten- 
tial. However, just as in Sec. 4.4 for the confining poten- 
tial, additional information is required to uniquely deter- 
mine the values of y^ax and ymm for any applied current 
/sD because a range of y^ax values exists for which corre- 
sponding values of y^m can be obtained that provide the 
correct quantum Hall voltage. 

5.2 Charge-Redistribution Potential and Confining 
Potentials 

What if the current exists along the device edges and 
within the interior? Many authors have considered this 
possibility [42, 57-65]. The case considered here in- 
volves a confining potential KCv) defined by Eqs. (6) 
and (7) and located on either side of the device, and an 
applied current-dependent charge-redistribution poten- 
tial Vr(y) defined by Eqs. (14), (15), (1 1), and (12). Note 
from Eqs. (15), (11), and (12) that Vr(y) is zero every- 
where if there is no applied current /sd because 
/c(rhs) = — /c(lhs) for that situation. 

Figure 8 is a schematic diagram of these confining 
and charge-redistribution potentials. The thick lines in- 
dicate the regions of occupied eigenstates, which extend 
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Fig. 8. Schematic drawing of the confining potentials VcCv) and the 
charge-redistribution potential V^iy). The thin lines are the two types 
of potentials. Thick lines indicate the regions of occupied eigenstates, 
which extend between yj^^ and y^nm- The origins of the confining 
potentials at ± A are greatly exaggerated for clarity, and j^ax and jmin 
are much farther from the device sides than in an actual example. 

between y^ax and ymm- The values of ymax and ymin are the 
same for the confining and charge-redistribution poten- 
tials. Increasing the current shifts the thick lines closer 
to one side of the device and farther away from the other 
side. The sign of the shift depends on the magnetic field 
and applied current directions. 

The electrical transport properties depend on the total 
potential Vt(j), but if we can unambiguously separate 
Vt(y) into the confining and charge-redistribution po- 
tential components then 



Vt(j) = Vc(j)+Vr(j). 



(18) 



Most of the information required for Eq. (18) is 
known. The potentials Vc(j) and Vr(y) are defined by 
Eqs. (6), (7), (11), (12), (14), and (15). For a given 
device we know the applied current /sd and the device 
width w. The current-independent parameters for the 
confining potential are zl = 0.5 |ULm and V^ = 0.75 V. If 
the geometry factor G of the charge-redistribution po- 
tential is assumed to be current-independent, then eval- 
uation of Eq. (16) with y^ax = 199.559 |ULm and w = 400 
|ULm, found for the case when /sd = |ulA in Sec. 4.3, 
gives the value G = 0.147. 

The problem remains, however, that just as in the 
cases of separate confining potentials and charge-redis- 
tribution potentials, there are still two free parameters 
jrnax and jmin- Ordinarily, it is impossible to uniquely 
determine the values of y^ax and y^^^ since the only other 
information is that the Hall voltage Vh is 



Vh = Rnhu = VtCjmin) - VtCjmax), 



(19) 



and there is a range of values for y^ax and y^in that satisfy 
this equation. It was possible, however, in an experi- 
ment of Cage and Lavine [20], to determine y^ax, and to 
thereby obtain a unique value of ymin from Eq. (19). 

Cage and Lavine [20] measured the quantized longi- 
tudinal voltage drops V, along a GaAs/AlGaAs device 
on the / = 2 plateau at 12.3 T for high currents in the 
breakdown regime [6, 7, 19, 66-70], and deduced the 
maximum electric field E^ax from a quasi-elastic inter- 
Landau level scattering model [18-20, 71]. The result 
was 

£max = 1.1 X 10' V/m at /sd = 215 |ulA (20) 

to excite the lowest V^ quantum voltage. 

Since E(y) = - VV(yX it follows from Eq. (18) that 



£t(jmax) = £c(jmax) + ^rCjmax)- 



(21) 



It is clear from Fig. 8 that E^ax will occur at the side of 
the device where Vdy) and Vr(y) combine to yield the 
largest value. This is at y^ax', hence 



^tCjmax) - ^m 



(22) 



We can therefore use the electric field Eqs. (6), (17), 
(21), and (22) to determine y^ax, and then the potential 
Eqs. (6), (7), (14), and (15), plus the quantum Hall 
voltage Eq. (19), to obtain ymin for the device of Ref. 
[20]. (Note that changing the values of y^ax and jmin also 
alters the values of /c(rhs) and /c(lhs), and thereby the 
value of /r.) There are now no free parameters, and one 
can obtain unique solutions to the total potential Vt(y) 
and to other transport properties. 

5.3 Results for the Total Potential 

The calculations for a 400 |ULm wide device on the 
i = 2 plateau at 12.3 T were done in Ref. [17]. The 
values of yrnax and ymm at 215 |ulA were 199.599 |ULm and 
— 199.515 |ULm, respectively, which are 55 and 66 cy- 
clotron radii away from the device edges. Therefore, 
there are no skipping orbits. Note from Sec. 4.3 that y max 
and jmin were both 60 cyclotron radii away from the 
edges when /sd = |ulA; thus y^ax increased by only 40 
nm between /sd = |ulA and 215 |ulA. 

Cage and Lavine [17] also did the calculations for 
/sd = 25 |ulA, which is a typical current in quantized Hall 
resistance standards measurements. They used a linear 
interpolation of y^ax between the values for /sd = |ulA. 
For convenience, their plot of the total potential Vfy) is 
reproduced in Fig. 9. This predicted potential is in ex- 
cellent agreement with the contactless experimental 
measurements shown in Fig. 6 of Fontein et al. [13], 
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Fig. 9. The total potential KCy) at /sd = 25 fjiA for the i-2 plateau. 
The conducting states of the potential extend between 
Jmin^- 199.554 |jLm and yniax = 199.564 jjim. The quantum Hall 
voltage is 0.323 V. 

which verifies that this model of a total potential com- 
posed of confining potentials plus a charge-redistribu- 
tion potential is reasonable. 

The model of Cage and La vine [17] appeared to not 
conserve charge. Slightly more electrons were redis- 
tributed towards the + y side of the device than were 
removed from the — y side. Therefore, there seemed to 
be an unaccounted excess of electrons. They pointed 
out that charge conservation could be accomplished by 
adjusting the origin to the right until the total potential 
VtCy) and the charge-redistribution function 



im^ 



,S^(3,)=_,:^^y^(3,) 



(23) 



were self-consistent. They did not add this complica- 
tion, however, since the potential distribution was al- 
ready in good agreement with the experiment of Fontein 
et al. [13]. The charge is actually conserved in their 
model without adjustment of the origin because eigen- 
states of the Landau levels in the 2DEG become occu- 
pied/unoccupied by electrons on each side of the device 
near y^ax and ymin rapidly tunneling from/to the ionized 
donor atoms in the AlGaAs layer located above the 
2DEG and maintaining charge equilibrium, even at low 
temperatures [72] ; thus the net charge within the GaAs 
and AlGaAs layers is always zero. 

Once VtCy) = K(y) + K(y) is known, the electric 
fields £cCy) = - VKCy) and E,{y) = - VVr(y) can be 
determined from Eqs. (6) and (17). The current density 
Jt(y) for electrons moving in the positive x direction is 
then 

Jfy) = a^Eiy) = ^ [E^(y) + Er(y)l (24) 



]m: 



I(y)= /t(j)dy 



Rh ' 



(25) 



where 



and 



/SD= /tCy)dy = /CyTnax)+/Cymin), 

yiT 



A/ = /(y2)-/(yi). 



(26) 



(27) 



Cage and Lavine [17] also determined the percentage 
of the total current in each of 20 equal segments across 
the 400 |ULm wide device. Their results for /sd = 25 |ulA 
and 215 |ulA for the / = 2 plateau at 12.3 T are repro- 
duced in Fig. 10. The current distributions are nearly 
symmetric across the device, and are virtually identical 
between 25 |ulA and 215 |jlA. 

The edge-channel current could only be along the 
right-hand side segment in Fig. 10 for these current and 
magnetic flux density directions, but 70 % of the cur- 
rent is in the 19 segments to the left of this segment. 
Also, the edge-channel current would by necessity be 
within the confining potential regions, but 97 % of the 
applied current is in the region between — A and A, 
which is outside of the confining potentials. Therefore, 
a large fraction of the applied current is in the device 
interior. (This conclusion also follows directly from the 
experimental potential distributions of Fontein et al. 
[13], whose measurements are independent of any 
model assumptions.) 
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Fig. 10. Percentages of the total current in twenty equally spaced 
segments across the device width for /sd = 25 |jlA and 215 |jlA. The 
distributions are identical for the two currents. 
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A plot for the current distribution at /sd = |ulA is not 
shown here, but there would be — 0.81 fxA and + 0.81 
|ulA in the left-hand and right-hand side segments, re- 
spectively, because /c(rhs) = — /c(lhs) = 0.81 |ulA. No 
current exists in the other eighteen segments for this 
case. 



6. Discussion and Conclusions 



current is in the region between — A and A where the 
edge-channel current cannot exist. Therefore, nearly all 
of the applied current is within the device interior, and 
not at the device edges. This conclusion is independent 
of our particular definition of the potential distribution 
since an identical current distribution can be obtained 
directly from the data of Fontein et al. [13]. 

Acknowledgments 



The potential distribution across the width of a quan- 
tum Hall effect device has been modeled as a combina- 
tion of confining potentials along the device sides and a 
charge-redistribution potential. The confining potential 
is due to surface charges and charge-depletion near the 
device edges. The charge-redistribution potential arises 
from the applied current. Normally, there would be in- 
sufficient information to uniquely define the total poten- 
tial, but by using breakdown data at large applied cur- 
rents [17, 20] an example was found where the potential 
could be obtained with no free parameters. The resulting 
potential distribution is in excellent agreement with ex- 
periment [13]. 

Using this potential distribution, I predict that skip- 
ping orbits do not occur in quantum Hall devices. In- 
deed, the nearest current-carrying occupied states are 55 
to 60 cyclotron radii away from the physical device 
edges for currents smaller than the critical current for 
the onset of breakdown. 

Edge-channel states are a useful mathematical con- 
struct, especially for very small applied currents. How- 
ever, edge-state models only predict that the applied 
current is proportional to the difference in chemical 
potential across the device, not how the current is dis- 
tributed across the device [29]. One must know the 
potential distribution to obtain the current distribution. 

Can edge-channel states be distinguished from bulk 
states? Yes, if we define which states are at the "edges" 
and which are in the "interior" region. The edge states 
must be entirely within the confining potential regions 
because that is where the current is located when a 
device is placed in a magnetic field and there is no 
applied current. Therefore, let us assume that the "edge" 
states are composed of current-carrying occupied states 
of the two confining potentials on either side of the 
device and of the charge-redistribution potential, and 
that these "edge" states are in the two regions between 
A and j^ax and between — A and jmin- For the examples 
in Sec. 5.3 at 25 fxA and 21 5 |ulA we find that only 4.9 % 
and 2.7 % of these two applied currents are located 
between A and j^ax and only — 1.9 % and 0.3 % are 
between — A and jmm- (The — 1.9 % value means that 
the current is in the — x direction over this region of the 
device at 25 |ulA.) Furthermore, 97 % of the applied 
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